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■ Abstract. We present the time-averaged characteristics of the Crab pulsar in the 0.75-30 MeV energy window using data 
from the imaging Compton Telescope COMPTEL aboard the Compton Gamma-Ray Observatory (CGRO) collected over its 9 
year mission. Exploiting the exceptionally long COMPTEL exposure on the Crab allowed us to derive significantly improved 
COMPTEL spectra for the Crab nebula and pulsar emissions, and for the first time to accurately determine at low-energy 7-rays 
the pulse profile as a function of energy. These timing data, showing the well-known main pulse and second pulse at a phase 
separation of ~ 0.4 with strong bridge emission, are studied together with data obtained at soft/hard X-ray energies from the 
ROSAT HRI, BeppoSAX LECS, MECS and PDS, at soft 7-rays from CGRO BATSE and at high-energy 7-rays from CGRO 
EGRET in order to obtain a coherent high-energy picture of the Crab pulsar from 0.1 keV up to 10 GeV. The morphology of the 
pulse profile of the Crab pulsar is continuously changing as a function of energy: the intensities of both the second pulse and the 
bridge emission increase relative to that of the first pulse for increasing energies up to ~ 1 MeV. Over the COMPTEL energy 
range above 1 MeV an abrupt morphology change happens: the first pulse becomes again dominant over the second pulse 
and the bridge emission loses significance such that the pulse profile above 30 MeV is similar to the one observed at optical 

r> ' wavelengths. A pulse-phase-resolved spectral analysis performed in 7 narrow phase slices consistently applied over the 0. 1 keV 

- 10 GeV energy interval shows that the pulsed emission can empirically be described with 3 distinct spectral components: 
i) a power-law emission component (1 keV - 5 GeV; photon index 2.022 ± 0.014), present in the phase intervals of the two 
pulses; ii) a curved spectral component required to describe soft 100 keV) excess emission present in the same pulse-phase 
intervals; iii) a broad curved spectral component reflecting the bridge emission from 0.1 keV to ~ 10 MeV. This broad spectral 
component extends in phase over the full pulse profile in an approximately triangular shape, peaking under the second pulse. 
Recent model calculations for a three-dimensional pulsar magnetosphere with outer magnetospheric gap acceleration by Cheng 
et al. (2000) appear at present most successful in explaining the above complex high-energy characteristics of the Crab pulsar. 

Key words, pulsars: individual: PSR B0531+21 - Stars: neutron - supernovae: individual: Crab nebula - Gamma rays: 
observations - X-rays: stars 



1. Introduction 

The Crab pulsar (PSR B053 1 +2 1 ) has been studied extensively 
over the entire electromagnetic spectrum with pulse profiles 
dominated by two pulses, separated ~ 0.4 in pulse phase and 
approximately aligned in absolute phase over all wavelengths. 

After the first detections of pulsed emission in the 
X-ray regime by Fritz et al. (1969 ; ~ 1-13 keV) and 
Bradt et al. (1969| ; 1.5-10 keV) , the first significant detection 



of pulsed soft 7-ray emission was reported by Kurfess (1971 



100-400 keV). A great boost forward was made by the X-ray 
instruments aboard the OSO-8 (Pravdo & Serlemitsos 1981; 2- 



50 keV, revealing spectral variations of the pulsed emission as a 
function of pulse-phase), HEAO-1 (Knight 1982L 1 8-200 keV; 
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introducing the possible existence of at least 2 pulsed emission 
components: one associated with the 2 main peaks and one with 
the bridge, the interva l between the two peaks) and HEAO- 
2 (Einstein) satellites ( friarnden & Seward 1984] , 0.1-4.5 keV; 
producing the first high-resolution (~ 4") image of the Crab 
nebula/pulsar in X-rays along with a high-res olution soft X-ray 
pulse profile). Recently, |Pravdo et al. (1997) presented the de- 
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tails of a pulse-phase-resolved spectral analysis of the pulsed 
emission in the 5-200 keV interval based on RXTE PCA and 
HEXTE data. They found systematic spectral changes in the 
photon power-law index as a function of pulse-phase across 
the interval of the pulsed emission. Their work is confirmed 
by the findings presented by Massaro et al. (2000), who used 



data from the narrow field instruments aboard BeppoS AX (0.1- 
300 keV). These authors also made an attempt to disentangle 
the two emission components, introduced by Knight (1982), 
assuming for one component (reflecting the emission in the 2 
peaks) the shape of the optical pulse profile and for the other 
component (the bridge emission) a shape from an analytical 
model with adjustable parameters. 

At medium 7-ray energies the first detection of pulsed 
radi ation (0.6-9 MeV) from the Cra b pulsar was reported 
by |Hillieretal. (1970"| se e also frValraven et al. ( 1 975) , 



Graser & Schonfelder (1982) Mahoney et al. (1984) and 



Agrinieret al. (1990) 

In the high-energy 7-ray domain (J> 30 MeV) the first in- 
dications for pulsed emission from the Crab pulsar were ob- 
tained from data collected by balloon-borne spark chambers or 



gas Cherenkov detector systems (see e.g. Browning et al. 1971 
Albats et al. 19721 iParlier et al. 19731 pVkBreen et al. 19731). A 



big step forward in this energy range was made by the SAS-2 
spark chamber experiment in the early seventies. In these data 
(20 MeV - 1 GeV) significant pulsed emission was reported 
by [Kniffen et al. (1974)| and [Thompson et al. (1977j The most 
detailed early information on the pulsed high-energy 7-ray 
properties of the Crab pulsar was, however, provided by the 
data from the European COS-B satellite ( |3ennett et al. 1977 



Wills et al. 1982| pearet al. 1987[ ). Significant bridge emis- 
sion was discovered in the combined COS-B Crab dataset, and 
the spectral characteristics of the pulsed and unpulsed (nebula) 
emission turned out to be quite diverse. Also, the alignment of 
the main pulse (PI) from the radio regime up to high-energy 
7-rays was shown by Wills et al. (1982). 

The launch of the Compton Gamma-Ray Observatory 
(CGRO; 20 keV-30 GeV) in April 1991 brought about an 
enormous improvement in the statistical quality of the 7- 
ray data. During its exceptionally long lifetime of more than 
9 years Crab pulsar data were collected by the Energetic 
Gamma-Ray Experiment (EGRET; 20 MeV - 30 GeV) show- 
ing clearly changing spectral behaviour as a function of 
pulse-phase (Nolan et al. 1993; Fierro 1995; Fierro et al. 
1998). The imaging Compton Telescope COMPTEL (0.75 - 
30 MeV) viewed the Crab each time simultaneously with 
EGRET. The results based on data from an early set of ob- 
servations performed during the first-year all-sky survey of 



the CGRO mission had been published by Much et al. (1995) 



and Carraminana et al. (1994). Qlmer et al. (1994) presented 
the first findings from the Oriented Scintillation Spectrometer 
Experiment OSSE (50 keV-10 MeV), and in a CGRO-paper 
(Ulmer et al. 1995) the first results from OSSE, COMPTEL 
and EGRET were combined. It was shown that the overall pulse 
phase-averaged spectrum is not well fitted by a single power- 
law, better with a broken power-law. Phase-resolved spectra 
were produed and fitted with broken power laws, selecting 
three phase intervals (peak 1, the bridge and peak 2). 



In this paper the final COMPTEL 0.75-30 MeV re- 
sults on the Crab pulsar/nebula are presented using i) data 
from all available COMPTEL Crab observations, ii) up- 
graded/improved response estimates and iii) improved data se- 
lection criteria. In order to obtain a broad high-energy pic- 
ture, we also consistently analysed in detail Crab data from 
the Italian/Dutch BeppoS AX satellite at lower energies (0.1- 
300 keV) and from EGRET at higher energies (30 MeV-10 
GeV). For some parts of the work we analysed additional data 
(e.g. soft 7-ray/X-ray data from the CGRO Burst and Transient 
Source Experiment, BATSE, and the ROSAT HRI; and for 
comparisons optical and UV data). 

2. Instrument description and observations 

COMPTEL is the imaging Compton Telescope aboard CGRO 
and operates in the 0.75-30 MeV energy range. Its detection 
principle relies on a two layer interaction: a Compton scatter in 
one of the 7 upper-detector (Dl) modules followed by a second 
interaction in one of the 14 lower-detector (D2) modules. Main 
measured quantities are the angles (x, tp) specifying the direc- 
tion of the scattered photon (from the interaction loci in Dl and 
D2) and the energy deposits in the D1/D2 modules where the 
interactions took place. From the last two quantities we can cal- 
culate the scatter angle Tp and the total energy deposit E'tot (for 
a full description see Schonfelder et al. 1993). Its energy reso- 
lution is 5-10% FWHM and, due to its large field of view of 
typically 1 steradian, it is possible to monitor a large part of the 
sky simultaneously with a position determination accuracy of 
~ 1°. The events are time tagged with an accuracy of 0. 125 ms 
and are converted to Coordinated Universal Time (UTC) with 
an absolute accuracy better than 100 /is using the on board mas- 
ter clock, serving also the other 3 CGRO instruments BATSE, 
OSSE and EGRET. 

In this study we selected all CGRO viewing periods for 
which the angle between the pointing axis (co-aligned with the 
COMPTEL/EGRET z-axis) and the Crab pulsar is less than 
30°. Details for each individual observation can be found in 
Table [j], which is self-explanatory. 

Because we have also included extensively archival 
EGRET data in the current study, a brief summary of this 
CGRO instrument is given as well. EGRET has a (gas-filled) 
sparkchamber and is sensitive to 7-rays with energies in the 
range 30 MeV to 30 GeV. In the mode used for most of the ob- 
servations the field of view is approximately 80° in diameter. 
Its effective area is approximately 1500 cm 2 between 200 and 
1000 MeV. The angular resolution is strongly energy depen- 
dent: the 67% confinement angles at 35 MeV, 500 MeV and 3 
GeV are 10?9, 1?9 and 0?5 respectively. The energy resolution 
AE/E is — 20% (FWHM) over the central part of the energy 
range. The relative timing accuracy is 8 /is and the absolute 
accuracy is better than 100 /is. For a detailed overview of the 
EGRET detection princ iple and instrument characteristics, see 



Thompson et al. (1993) 



The X-ray data most extensively used in this work had 
been collected with the 4 narrow field instruments aboard 
BeppoSAX: the low-energy (0.1-10 keV) and medium energy 
(1.6-10 keV) concentrator spectrometers, LECS and MECS 
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Table 1. COMPTEL observation summary with PSR 
B0531+21 

within 30° from the pointing axis 
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Start Date 
TJD f 



End Date 
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respectively, the High-Pressure Gas Scintillation Proportional 
Counter, HPGSPC (4-60 keV) and the Phoswich Detector 
System PDS sensitive in the 15-300 keV energy range. Detailed 
instrument descript ions for the 4 narr o w field instruments 
can be found in |Parmar et al. (1991% |Boella et al. (1997), 
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Fig. 1. Pulse profiles of PSR B0531+21 as observed by CGRO 
COMPTEL in the 0.75-30 MeV energy interval (top) and 
CGRO EGRET (bottom) for energies > 10 MeV. The bound- 
aries of the pulse phase intervals defined by Fierro et al. (1998) 



are indicated by the dotted vertical lines. Notice the morphol- 
ogy change: at energies between 0.75-30 MeV the second peak 
dominates, while at energies above 100 MeV the first peak 
dominates. Moreover, in the 0.75-30 MeV interval there is con- 
siderable bridge emission in between the two peaks, which is 
hardly present at energies above 100 MeV. 



Manzo et al. (1997)| and [Frontera et al. (1997)| for the LECS, 
MECS, HPGSPC and PDS respectively. 



3. Timing analysis 

The first step in the timing analysis is to subject the events reg- 
istered during an observation to an event selection filter. In the 
case of COMPTEL the most important selection parameters are 
the t ime-of-flight TOF, the p ulse shape discriminator PSD (see 
e.g. |Sch6nfelder et al. 1993 1, the "spatial" parameters x,^,Tp 
and the total energy deposit E tot . Given the a priori known po- 
sition of the Crab pulsar it is possible to determine for each 
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Table 2. Phase component definitions for the Crab pulsar 
adopted in this study (see also Fig. [l]) 



Component 


Abbreviation 


Phase interval 


Width 


Leading Wing 1 


LW1 


0.88- 


0.94 


0.06 


Peak 1 


PI 


0.94- 


1.04 


0.10 


Trailing Wing 1 


TW1 


0.04- 


0.14 


0.10 


Bridge 


Bridge 


0.14- 


0.25 


0.11 


Leading Wing 2 


LW2 


0.25- 


0.32 


0.07 


Peak 2 


P2 


0.32- 


0.43 


0.11 


Trailing Wing 2 


TW2 


0.43- 


0.52 


0.09 


Off Pulse 


OP 


0.52- 


0.88 


0.36 


Total Pulse 


TP 


0.88- 


1.52 


0.64 



event the so-called (p geo angle, i.e the angle between the scat- 
tered photon and the source. The Tp angle provides an equiva- 
lent measure of this angle, but now only based on the energy 
deposits in both detector layers. The difference angle ip arm = 
Tp — ip geo is called the Angular Resolution Measure (ARM) and 
forms the base of the spatial response of COMPTEL and its 
distribution is narrowly peaked near tp Mm = with asymmetric 
(energy dependent) wings. The definite and significant timing 
signature of the Cr ab pulsar in the C OMPTEL energy range 
(Much et al. 1995; Much et al. 1997 ) provides a very usefull 
tool to determine the optimum event parameter windows for 
celestial sources. In this study we have determined and used 
subsequently the optimum (total energy deposit dependent) pa- 
rameter windows for the TOF and PSD. The optimum windows 
for the (p m - m angle turn out to be asymmetric around and a 
function of total energy deposit (as expected). Finally, for the 
given combination of viewing periods (see Table [j]) we com- 
pared the measured Tp distribution, dominated by background 
photons, with the distribution expected for a point source at the 
Crab position. This allows for a determination of the optimum 
window for selection on Tp (a function of the total energy de- 
posit). 

Events, not vetoed by any of the 4 anti-coincidence domes 
surrounding the COMPTEL detector layers, and having a very 
low probability of originating from the Earth ("Earth Horizon 
angle selection") and passing through our optimized event se- 
lection windows are finally used in the timing analysis. 

The event arrival times (at the spacecraft) of these accepted 
events are converted to arrival times at the barycentre of the so- 
lar system (SSB) using the JPL DE200 solar system ephemeris 
and the Crab pulsar position. This process yields SSB arrival 
times with an absolute timing accuracy of less than 100/^s. 
Subsequent phase folding using the Crab pulsar ephemeris data 
(CGRO timing database; Arzoumanian 1992) yields the pulse 
profile. Because the Crab pulsar shows a lot of timing noise 
(young pulsar) the ephemerides have only a limited validity in- 
terval. Therefore, in the phase folding process we used different 
timing solutions (ephemerides) for (almost) each observation 
given in Table |l[ 

The derived final - combining all observations given in 
Table 0- COMPTEL pulse profile over the total energy range 
0.75-30 MeV is shown in Fig. [j]. The two expected peaks, sep- 
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Fig. 2. Pulse profiles (double cycles) of PSR B0531+21 as 
measured by CGRO COMPTEL in 4 differential energy inter- 
vals: 0.75-1, 1-3, 3-10 and 10-30 MeV. Typical error bars are 
shown near phase 0.6. A clear morphology change of the pulse 
profiles is visible: below 10 MeV the second peak (near phase 
0.4) dominates, and emission in the "bridge" phase interval is 
significant, while above 10 MeV the first peak (near phase 0.0) 
dominates with strongly reduced "bridge" emission. 



arated ~ 0.4 in phase, with intense emission between the peaks 
(bridge emission) are visible with high statistics. The total ex- 
posure has increased by a factor ~ 5 compared to the last pub- 
lished total COMPTEL profile of the Crab (Much et al. 1995). 
For comparison we also show in Fig. |l] the CGRO EGRET 
profile for energies above 100 MeV for which we analyzed 
archival Cycle to VI EGRET viewing periods with the Crab 
pulsar within 35° from the pointing axis, and in which the spark 
chamber was switched on. The same ephemerides have been 
used as for the COMPTEL data. The differences in morphol- 
ogy between the COMPTEL and EGRET profiles are evident. 
Following Fierro (1995) we used the EGRET profile shape to 
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select narrow phase intervals for phase resolved spectroscopy 
studies (see below). The intervals are shown in the figure and 
given in Table ^. Fig. ^| shows the pulse profiles (double cy- 
cles for clarity) in the 4 "standard" COMPTEL energy intervals 
(0.75-1, 1-3 3-10 and 10-30 Me V ). The significances applying 
the Z^-test (Buccheri et al. 1983) with 8 harmonics on the un- 
binned set of pulse-phases are (expressed in Gaussian sigma's) 
20.0cr, 31.7(7, 18.6cr and 10.9(7, respectively. Comparing the 
profiles with those presented in Much et al. (1997) shows the 
enormous increase in statistical quality, especially for energies 
above 3 MeV, due to the longer exposure in combination with 
our improved event selection procedures. We see the morphol- 
ogy change from the COMPTEL to the EGRET profile in Fig. 
occurring over the COMPTEL energy range: below ~ 10 MeV 
the second peak (near phase 0.4) dominates the first peak (near 
0), and the bridge emission is intense, while above ^10 MeV 
the first peak dominates and the bridge emission is strongly re- 
duced. 

The 7-ray profiles in Fig. [j] and Fig. || are time-averaged 
profiles, compiled over many years. Before analysing the pro- 
files further, we first verified the long-term stability of the 7-ray 
signature (flux and pulse shape). 

4. Long-term variability 

4. 1. Flux variations: total pulsed flux in the 1-10 
MeV range 

We studied the time variability of the emission from the Crab 
pulsar by determining the "pulsed" flux in differential energy 
windows as a function of time. The "pulsed" flux has been de- 
rived from the number of excess counts in the Total Pulse phase 
interval (see Table |^) on top of the average emission level in the 
Off Pulse phase interval. For the latter interval we assume that 
the emission originates from the nebula only, although a DC- 
contribution from the pulsar can not be ruled out. We show as 
examples the 1-3 and 3-10 MeV results. These results over the 
more than 9 year baseline (April 1991/ May 2000) are shown 
in Fig. |[ For the 1-3 MeV energy interval typical integration 
times are 2/3 weeks, while in the 3-10 MeV interval longer in- 
tegration times are used given the strongly reduced statistics. 
The xt values for the fits assuming a constant flux level are 
typically ~ 0.5 for both energy ra nges, indicatin g that there 
is no evidence for time variability. Fierro (1995) also studied 
the Crab long-term variability for energies above 100 MeV and 
also concluded that the emission from the Crab pulsar is stable. 

4.2. Pulse shape variations: P2/P1 ratio in the 1-3 
MeV range 

The time variability of the pulse shape was investigated in the 
1-3 MeV (best statistics) energy window by determining the 
P2/P1 ratio for each observation. This ratio is derived by mea- 
suring the number of excess counts in the P2 and PI phase 
intervals (see Table ^ on top of the level in the Off Pulse in- 
terval. This ratio is shown in Fig. The xl value for a fit 
assuming a constant P2/P1 ratio is ~ 0.8. Therefore, there 
is no indication for a time dependency of the pulse shape 
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Fig. 3. "Total Pulsed" flux from the Crab pulsar as a function 
of time in the 1-3 MeV (top) and 3-10 MeV (bottom) energy 
intervals. The ±1(7 uncertainty intervals assuming a constant 
flux are indicated by the shaded regions in both figures. The x1~ 
values for the fits assuming the flux being constant are ~ 0.5, 
thus there is no indication for "Total Pulsed" flux variability in 
both the 1-3 and 3-10 MeV energy intervals. 
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Fig. 4. The P2/P1 flux ratio of the Crab pulsar in the 1-3 MeV 
energy interval as a function of time. The ±1(7 uncertainty in- 
terval assuming a constant ratio is indicated by the shaded re- 
gion. The P2/P1 ratio is consistent with being constant. 



at medium energy 7-ray energies consistent with the findings 
prese nted by jCarraminana et al. (1994)[ For energies above 30 
MeV |Fierro (1995f studied the long-term temporal variation 
of the P2/P1 ratio and found also no evidence for a (sys- 
tematic) variation over the Cycle 0-III EGRET observations. 
Tompkins et al. (1997) came to a similar conclusion using an 
extended EGRET data base including also Cycle IV- V. 
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5. Pulse profiles of PSR B0531+21 from 0.1 keV 
up to 10 GeV 

The 7-ray pulse profiles in Fig. || show that the pulse mor- 
phology changes significantly over the COMPTEL energy win- 
dow (0.75-30 MeV), i.e. the emission spectra vary significantly 
with phase. Phase-resolved spectral analyses have earlier been 
performed at X-ray and 7-ray energies for different data sets 
and/or different narrow energy intervals. However, for each 
study, different phase selections have been made such that a 
consistent full high-energy picture of the Crab pulsar can not be 
compiled from published results. Therefore, we extended our 
energy window by analysing consistently not only the CGRO 
EGRET (30 MeV - 10 GeV) high-energy 7-ray data, but also 
X-ray/soft 7-ray data from the ROSAT HRI (0.1-2.4 keV), 
BeppoSAXLECS (0.1-10 keV), MECS (1.6-10 keV) and PDS 
(15-300 keV), and CGRO BATSE (20 keV- 1 MeV). 

The (on board folded) data from CGRO BATSE overlap in 
energy with the data from CGRO OSSE for which results have 
already been published by Ulmer et al. (1994,1995). However, 
due to the enormous exposure in the co-added BATSE data the 
statistics are much better than can be obtained in the combined 
OSSE Crab observations. Especially above ~ 220 keV (e.g. 
Ulmer et al. 1994, Fig. 2) where the OSSE data have low sta- 
tistical quality, the BATSE profiles are superior. For the high 
statistics OSSE data (below ~ 220 keV) we verified that the 
profiles are consistent in shape with those obtained by us using 
BeppoSAX PDS and CGRO BATSE data. 

In the next subsections details are given about the compila- 
tion of Crab pulse profiles over the energy range 0. 1 keV to 10 
GeV. 

5. 1 . ROSAT HRI 0. 1 -2. 4 keV pulse profile 

The soft X-ray ROSAT HRI data were collected during an ob- 
servation of the Crab pulsar/nebula performed from 4 March 
1995 to 15 March 1995 yielding a net exposure time of 
7.98 ks (HEASARC Online Service; observation identifier 
RH400639N00). Because the data are spread over 1 15 different 
orbital intervals over the 1 1 day observation period the consid- 
erable ROSAT clock drift will result in a messy pattern when 
combining the pulse phases from the entire observation. We 
could identify 4 consecutive sets of orbital intervals in which 
the observed pulse profile is stable. The 0. 1 -2.4 keV pulse pro- 
file shown in Fig. gp was obtained cross-correlating 3 of the 4 
profiles with the profile chosen as template, correcting for the 
observed mutual phase shifts and fixing the zero phase at the 
centre of the main peak. 

5.2. BeppoSAX 2.4-1 00 keV pulse profiles 

The BeppoSAX LECS, MECS and PDS data have been 
collected during a calibration observation of the Narrow 
Field instruments aboard BeppoSAX performed on 25-26 
September 1999 yielding (screened) effective exposure times 
of 7.75 ks, 32.6 ks and 30.7 ks for the LECS, MECS (unit- 
2) and PDS clusters A & B, respectively (data retrieved 
from archive maintained by BeppoSAX ASI Science Data 
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Fig. 5. High-energy pulse profiles of PSR B0531+21 from 0.1 
keV up to 10 GeV. Data have been used from the following in- 
struments: a) ROSAT HRI (0.1-2.4keV), b) BeppoSAX MECS 
(2.4-10 keV) , c) BeppoSAX PDS (20-100 keV), d,e) CGRO 
BATSE (100-3 15 keV & 3 15-750 keV), f,g) CGRO COMPTEL 
(0.75-10 MeV & 10-30 MeV) and h) CGRO EGRET (> 30 
MeV). The morphology change of the profiles as a function of 
energy is striking. 



Center at http://www.asdc.asi.it/bepposax/; Observation Codes 
20795007 & 207950071). In Fig. |5jb and c the Crab pulse pro- 
files are shown as observed by the MECS in the 2.4-10 keV 
energy window and by the PDS in the 20-100 keV energy win- 
dow, respectively. 

5.3. CGRO BATSE 100-745 keV pulse profiles 

In the hard X-ray/soft 7-ray band (0.05-1 MeV) we have used 
archival data from the CGRO BATSE Large Area Detectors 
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collected during observations performed between MJD 48392 
and 50273 in the onboard folding mode (CGRO Archive main- 
tained by HEASARC at ftp://cossc.gsfc.nasa.gov/compton/data 
directory batse/pulsar/onboardjolded/crab/). Typical integra- 
tion times were 2/3 weeks per included observation (64 obser- 
vations have been used in this study). The profiles had been 
produced in 64 bins per cycle in 16 different energy channels 
for each individual observation run. We determined the shifts 
of the pulse profiles of the individual observation runs with re- 
spect to the profile obtained during the observation run 48392- 
48406 in channel 9 (~ 165-230 keV) by cross-correlation. 
Applying the shifts in the combination of the pulse profiles and 
putting the first peak at phase yields high quality pulse pro- 
files in the 20 keV - 2 MeV range. In Fig. and e the profiles 
are shown for the 100-315 keV (channels 7-10) and 315-745 
keV (channels 11-13) energy windows, respectively. 

5.4. CGRO COMPTEUEGRET pulse profiles 

At medium energy 7-rays the CGRO COMPTEL pulse pro- 
files derived in this work are shown for the 0.75-10 and 10-30 
MeV energy windows in Fig. |5]f and g, respectively. Note that 
the COMPTEL profiles have a large non-zero offset. Finally, in 
Fig. |h the CGRO EGRET pulse profile is given for energies 
above 30 MeV (we used data from Cycle 0-VI observations, 
retrieved from the CGRO Archive maintained by HEAS AR C; 



ftp://cossc.gsfc.nasa.gov/compton/data/egret/highJevelfy . 



6. P2/P1 and Bridge/P1 ratios as a function of 
energy 

From the (high-energy) pulse profile compilation shown in Fig. 
H we can immediately observe some striking features. The sec- 
ond peak (near phase 0.4) becomes more and more pronounced 
for increasing energies. However, above ~ 10 MeV the first 
peak becomes dominant again. The "Bridge" emission seems 
to show a similar behaviour as the second peak. In a more quan- 
titative evaluation of this morphology change of the profile as a 
function of energy we determined the intensity ratios for P2/P1 
and Bridge/Pi as a function of energy over the entire range 0. 1 
keV to 10 GeV, adopting the phase interval definitions of Table 
||. The pulsed emission in each interval has been separated from 
the underlying nebula/DC emission by subtracting the (prop- 
erly scaled) emission from the OP phase interval. The results 
are visualized in Figs. || & for the P2/P1 and Bridge/Pi ra- 
tios, respectively. In these plots we have also included the ratios 
derived from the optical profile in the 3800 to 6500 A wave- 
length interval obtained by Much et al. (2000) using the UCL 
MIC detector as well as those in the far-ultraviolet (1 140-1720 
A) and near-ultraviolet (1600-3200 A), obtained from (time- 
tagged) data taken by the HST STIS instrument (Sollerman et 
al. 2000, Gull et al. 1998). In all optical ranges we again ap- 
plied the consistent phase interval definitions. 

The P2/P1 ratio as a function of energy (Fig. ^) gradu- 
ally increases from the optical wavelength range to ~ 1 MeV, 
followed by a rapid decrease in the 1-30 MeV interval (the 
COMPTEL energy window) towards a more or less constant 
value of ~ 0.5 for energies above 30 MeV (the EGRET en- 
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Fig. 6. P2/P1 ratio as a function of energy from optical wave- 
lengths up to high-energy 7-rays. Data from the following in- 
struments have been used: optical wavelengths, UCL MIC de- 
tector (star symbol), NUV/FUV HST STIS (star symbols); X- 
ray energies, ROSAT HRI (open square), BSAX LECS (filled 
square), BSAX MECS (open triangle); Hard X-rays/soft 7- 
rays, BSAX PDS (open circles), CGRO BATSE (filled cir- 
cles); Medium/hard 7-ray energies, CGRO COMPTEL (filled 
upwards pointing triangle), CGRO EGRET (filled downwards 
pointing triangle). The gradual increase of the P2/P1 ratio up to 
~ 1 MeV is striking, a sharp decline in the 1-30 MeV energy 
range follows and a recovery to the optical ratio value settles 
above ~ 30 MeV. 

ergy window). The Bridge/Pi ratio vs. energy (Fig. (7J) exhibits 
a very similar shape as for the P2/P1 ratio. However, the values 
of the latter ratio in the optical and high-energy 7-ray domains, 
become very small (0.017) indicating that the Bridge emission 
practically vanishes. It is only substantial in the ~ 1 keV to 
~ 10 MeV energy window in contrast with the emissions from 
the 2 peaks which are always present. This behaviour suggests 
that we are dealing with an emission component distinct from 
the emission from both peaks. This hypothesis (see e.g. also 



Knight 19821 Hasinger 1984,1985; Massaro et al. 2000) is fur- 
ther strengthened in the phase-resolved spectral analysis pre- 
sented in the next section. 

Similar analyses have been presented in the past by other 
authors (e.g. Toor & Seward 1977; Hasinger 1984,1985; Ulmer 
et al. 1994; Mineo et al. 1997; Massaro et al. 1998,2000), gen- 
erally over more restricted energy windows with poorer data 
coverage, and/or often using data of inferior statistical qual- 
ity. Particularly the "transition" region of the COMPTEL MeV 
window is now well covered for the first time. 

7. Spectral analysis 

In this section we will first present the spectra of the nebula 
emission and the Total Pulse emission (excess emission in Total 
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Table 3. COMPTEL spectra of the Crab nebula and pulsar 
(Total Pulse). Fluxes with la statistical uncertainties 
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Fig. 7. Bridge/Pi ratio as a function of energy from optical 
wavelengths up to high-energy 7-rays. See the caption of Fig. |] 
for the meaning of the symbols. The Bridge/Pi ratio is almost 
at optical wavelengths, but gradually reaches a maximum near 
1 MeV, followed by a drastic break in the 1-30 MeV energy 
range. Above ~ 30 MeV the Bridge/Pi ratio approaches the 
optical value of ~ 0.017. 



Pulse interval; see Table ||). Then we will show the results from 
the phase-resolved spectral analysis for the narrow phase in- 
tervals. As for the timing analysis, we did not limit ourselves 
to the analysis of the COMPTEL data, but we collected data 
over a very wide spectral band to derive a consistent overall 
high-energy picture of the Crab. Combining spectra derived by 
different instruments, we had to assess possible systematic ef- 
fects in the flux estimates and their impact on our analysis and 
conclusions. For this purpose, we used the nebula spectra to 
estimate possible inconsistencies. 

7. 1. Crab nebula spectrum 

To determine the Crab nebula spectrum from our COMPTEL 
(0.75-30 MeV) Cycle 0-IX observations (see Table [l]) we se- 
lected events recorded in the Off Pulse phase interval (see Table 
|J), assuming that any pulsar DC emission is negligible, and 
applied a maximum likelihood method using the spatial sig- 
nature of a point source in the (7p, ip aim ) plane as a function 
of measured energy. This work yielded an improved spectrum 
for the nebula emission in the COMPTEL energy range com- 
pared to the COMPTEL Crab nebula spectrum published ear- 
lier by van der Meulen et al. (1998). In the latter work a smaller 
database was used, as well as preliminary response characteris- 
tics, which have since been improved upon. The newly derived 
COMPTEL nebula spectrum is given in Table ||. A power-law 
fit to the COMPTEL nebula flux points results in a photon in- 
dex of 2.227 ± 0.013. In the same table also the statistical un- 
certainties (la) on the flux measurements are provided. In this 
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context we note that COMPTEL has an overall (systematic) 
uncertainty on its flux estimates of the order of 10 — 20%. 

To cover the neighbouring soft X-ray to soft 7-ray band, 
we derived the BeppoSAX LECS, MECS, HPGSPC and PDS 
nebula spectra applying the most recent (December 1999 is- 
sue for the LECS and HPGSPC and November 1998 issue for 
the MECS and PDS) response characteristics (sensitive area, 
energy redistribution matrices and spatial response). We fitted 
the Sept. 1999 Off Pulse Crab data from these four BeppoSAX 
instruments simultaneously over the full 0.1-300 keV energy 
range with an absorbed power-law model taking into account 
the mutual uncertainties in absolute flux calibrations by in- 
cluding in the fit three free relative normalization scale fac- 
tors (MECS scale factor fixed to 1). The energy of the LECS 
events used in the fit was constrained to the 0.1-4 keV win- 
dow. The best photon index and Hydrogen column density 
Nh, assuming solar abundances for the other elements absorb- 
ing soft X-rays, resulting from this fit are 2.145 ± 0.001 and 
3.61(2) x 10 21 cm~ 2 , respectively. In Appendix^ the system- 
atic uncertainties in the derived spectral characteristics in the 
X-ray/soft 7-ray band are discussed in detail comparing flux 
estimates obtained by different high-energy instruments. 

For the adjacent EGRET high-energy (30-30000 MeV) 7- 
ray range we derived (phase-resolved) Crab spectra using all 
available Cycle 0-IV archival EGRET data for which reli- 
able sensitivity estimates were available to us[] The method 
is equivalent to the spatial maximum likelihood analysis per- 
formed by Fierro (1995). The only difference with the latter 
work is that we now added Cycle-IV data. The EGRET spect ral 
data are claimed to be 10% accurate (Thompson et al. 1995). 



Due to the problem of gas aging, the spark chamber efficiency 
degraded significantly during the later Cycles of the EGRET CGRO 
observations. Recently, the EGRET team did present the energy and 
time dependent correction factors for the later Cycles (Gamma 2001 
Symposium, Baltimore). 
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Fig. 8. The Crab nebula spectrum from 
soft X-rays up to TeV 7-rays. The TeV 
data point near 1.6 x 10 5 MeV is taken 
from Oser et al. (2001) and the hatched 
bands represent the flux measurements 
and corresponding la uncertainty esti- 
mates at TeV energies (for references, 
see text). 
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Fig. 9. The Total Pulse emission of the 
Crab pulsar from optical wavelengths 
up to high-energy 7-rays. The nebula 
emission has been subtracted. The opti- 
cal spectral data (10~ 6 -10~ 5 MeV) are 
taken from Sollerman et al. (2000) and 
the TeV data point near 1.6 x 10 5 MeV 
from Oser et al. (2001). 



Energy [MeV] 



The Crab nebula spectrum from 0. 1 keV up to 50 TeV is 
shown in an E 2 x F representation in Fig. [| Included are the 
BeppoSAX, COMPTEL and EGRET spectra derived in this 
work, soft 7-ray spectral information from GRIS (0.02-1 MeV; 
Bartlett et al. 1994a) together with ground-based TeV data 
(STACEE-32 > 0.19 TeV, bseretal. 20011 HEGRA 1-20 TeV, 



Aharonian et al. 2000 Whipple 0.5-8 TeV, 


Hillas et al. 1998 


CANGAROO 7-50 TeV, Tanimori et al. 1998 


). Similarly to the 



spectrum shown in van der Meulen et al. (1998), but now more 
pronounced, we see a continuous and smooth decrease from 
soft X-rays up to medium energy 7-rays, irrespective of the un- 
certainties in the absolute sensitivities of the instruments, fol- 
lowed by a steep gradient beyond ~ 30 MeV to ~ 300 MeV. 
Above ~ 300 MeV an additional emission feature seems to 
emerge reaching a maximum between 10 and 100 GeV, a win- 



dow which is not yet accessible for space-borne and ground- 
based experiments. For an interpretation of this spectral shape, 
see e.g. de Jager et al. (1996) 



7.2. Crab Total Pulse spectrum 

In the 0.75-30 MeV energy range we determined the pulsed 
flux values from all CGRO COMPTEL observations (see Table 
|l]) in two distinct manners. In the first method, the number of 
(pulsed) excess counts in the broad Total Pulse interval (see 
Table |J) above the mean level in the Off Pulse interval is deter- 
mined as a function of measured energy. These excess counts 
are then converted to flux values using the COMPTEL re- 
sponse and the total exposure. The second approach is based 
on the maximum likelihood method in the spatial domain as 
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introduced in Sect. 7.1 . Applying the latter approach for pulse- 
phase selected events and subtracting the properly scaled Off 
Pulse contribution (containing only the DC/nebula source with 
a point source signature) yields the pulsed fluxes as a function 
of energy for the selected pulse-phase intervals. We verified for 
the Total Pulse interval, having the best statistics, that the fluxes 
derived from both methods are compatible within 5 — 20%, giv- 
ing a measure of the systematic uncertainties. The COMPTEL 
Total Pulse fluxes from the spatial analysis (remember, nebula 
emission subtracted) are also included in Table ^ and shown in 

For comparison, we included in Fig. |9j our derived Total 
Pulse Crab EGRET spectrum (30 MeV- 10 GeV, Cycle 0-IV) 
and the published Total Pulse spectra from GRIS (20 keV - 1 
MeV: Bartlett 1994a) and OSSE (50 keV - 0.59 MeV: Ulmer et 
al. 1994). In the latter two publications slightly different phase 
intervals have been used to derive the Total Pulse spectrum. 

In the COMPTEL Total Pulse spectrum a feature becomes 
apparent: the high flux value in the 10-15 MeV interval, con- 
sistently derived in both the timing and spatial methods. A 
response anomaly is excluded, e.g. the nebula 0.75-30 MeV 
spectrum (cf. Fig. ||) exhibits a very smooth behaviour over 
its entire range, nor has such an effect been seen in other 
COMPTEL analyses. Fitting the COMPTEL flux points with 
a power-law spectral model, excluding the deviant 10-15 MeV 
flux point, yields a good fit with photon index of 2.35 ± 0.06 
(xt — 0-50 for 9 d.o.f.), connecting smoothly at both ends to 
the GRIS/OSSE and EGRET flux measurements. The excess 
flux in the 10-15 MeV interval above this power-law model 
fit reaches a significance of 3.5<r. Including the deviant 10-15 
MeV flux point in the spectral fit yields a worse fit (xt = 1.42 
for 10 d.o.f) with a power-law index of 2.24 ±0.04, which does 
not connect smoothly to the neighbouring measurements, par- 
ticularly to EGRET. The 10-15 MeV excess flux above this fit 
has a significance of only 2. 0<r. We do therefore not regard the 
10-15 MeV flux enhancement as a firm detection of a new spec- 
tral feature, and have no possible astrophysical interpretation, 
but we find it interesting to note that contributions to this flux 
enhancement appear to come from those (narrow) phase inter- 
vals in which a broad spectral (Bridge) component dominates 
the spectrum (see Sect. 7.3). If genuine, it could therefore be 
related to this spectral component. 

For the BeppoSAX LECS, MECS and PDS we determined 
the number of (pulsed) excess counts above the mean level in 
the Off Pulse phase interval, similarly to the first COMPTEL 
method. These excess counts have been converted to flux mea- 
sures applying the most recent response characteristics assum- 
ing intrinsic power-law type emission absorbed in a column of 



density N H = 3.61(2) x 10 21 cm' 2 (see Sect. 7.1). In this way 



we obtained the Total Pulse spectrum over the 0.1-300 keV en- 
ergy interval, which is also shown in Fig. 

We augmented the energy coverage by including the 
pulsed spectra obtained at optical/NUV/FUV wavelengths by 
Sollerman et al. (2000). In this E 2 x F representation, the 
total pulsed emission shows a gradual increase from the opti- 
cal range towards a plateau of maximum luminosity extending 
from ~ 10 keV to ~ 1 MeV. Beyond ~ 1 MeV the emission 
softens until ~ 70 MeV, above which a second plateau appears 



with an emission spectrum having a photon power-law index 
close to 2. Between 4 and 10 GeV the pulsar spectrum appears 
to break/soften drastically to account for the recently derived 
upper limits for pulsed emission at TeV 7-ray energies (see e.g. 

Borione et al. 1997[ A.haronian et al. 1999 



Vacanti et al. 1991 



Lessard et al. 200C and Dser et al. 2001 ). Its spectral behaviour 
is completely different from that of the nebula (cf. Fig. |]). 
Notice the dominance of the nebula emission component over 
the pulsed emission component for energies below ~ 100 
MeV and above ~ 10 GeV, comparing Figs. |8| & ^. Only in a 
small window at high-energy 7-rays between ~ 100 MeV and 
~ 10 GeV the pulsed component exceeds the underlying neb- 
ula component. This Total Pulse spectrum is clearly complex. 
For detailed theoretical interpretations it is important to disen- 
tangle first this total spectrum in contributions from different 
phase components of the pulse profile. 



7.3. Spectral behaviour in the narrow pulse-phase 
intervals 

A similar broad-b and spectral analysis (0.1 keV - 10 GeV), as 
presented in Sect. 1J1 for the Total Pulse interval, has been per- 
formed for the 7 narrow pulse-phase intervals defined within 
the Total Pulse phase interval (see Table |J). We analyzed 
data from the BeppoSAX LECS, MECS and PDS and CGRO 
COMPTEL and EGRET instruments applying identical phase 
window selections. Since we planned to make empirical fits 
to the multi-instrument spectra, we wished to avoid being too 
sensitive to the systematic discrete jumps in the overlapping 
spectra of the BeppoSAX LECS, MECS and PDS and GRIS, 
as shown in Fig. A. 1 . In the spectra and analysis presented here, 
we added a representative 10% systematic flux uncertainty to 
the statistical uncertainty in each flux measurement involved 
(for COMPTEL 15% was used). We also repeated the total 
analysis normalizing the BeppoSAX data on the much lower 
GRIS value, as well as on the average normalization value of 
the BeppoSAX instruments. Our results are not sensitive to 
these different normalizations. 

The spectr al re sults for the seven different intervals are 
shown in Fig. 1(J. This compilation clearly shows that the 
spectral shape varies strongly with pulsar phase. For example, 
the spectra of the PI and the Bridge intervals differ dramat- 
ically. Emission from the latter is hardly discernible at ener- 
gies below 1 keV and above 100 MeV; this emission is con- 
fined roughly between these energies in a broad "bump" shaped 
fashion in this E 2 x F representation. On the other hand, the 
emission in the PI interval remains very strong in the 7-ray do- 
main above 1 MeV (COMPTEL and EGRET data), exhibiting a 
power-law photon distribution up to a spectral break at GeV en- 
ergies. Extrapolation of such a power-law spectral shape to X- 
ray energies, reveals a narrower "bump" shaped excess above 
the power-law extrapolation, with a maximum power output 
well below 50 keV. The P2 spectrum is rather similar to the 
spectrum of PI, but a spectral component similar in shape 
to that of the Bridge interval seems to enhance the P2 spec- 



~ The phase-resolved spectral data can be retrieved from 
\ittp://wsl 3. sron.nl: 8080/personal/kuiper/data 
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Fig. 10. The high-energy emission of the Crab pulsar in the 7 narrow pulse-phase intervals (see Table Q) from 0. 1 keV up to 
10 GeV. Two spectra (for the TW1 and LW2 phase intervals) are displayed twice, to facilitate a better visual comparison of the 
different spectra (see discussion). See Fig. ^ for the meaning of the symbols of the data points. The different contributions (PI 
modified power-law model, dash-dot-dot-dot line; Bridge modified power-law model, dash-dot line; power-law model, dash line) 
to the composite model fits (dotted lines) have all been superimposed for the best-fit scaling parameters (histogram) shown in 
Fig-0 



trum at MeV energies, relative to the PI spectrum. Finally, the 
Trailing Wing 1 (TW1) and Leading Wing 2 (LW2) spectra 
are amazingly similar in shape, and appear to be some mix- 
ture of shapes of PI and the Bridge intervals. This should not 
be a surprise, because the adopted separation in phase intervals 
(Table will most likely not coincide exactly with genuine 
physical components (different dominating production mech- 
anisms and/or production sites in the Crab magnetosphere). 
However, the vastly different spectral behaviour exhibited in 
the Bridge phase interval suggests a physically distinct emis- 



sion component as proposed earlier by e.g. Knight (1982) and 
ftasinger(1985j[ 



7.4. Parametrization of the emission in the narrow 
pulse-phase intervals 

Exploiting our high statistics and eight-decades wide high- 
energy phase-resolved spectra, we made an attempt to empiri- 
cally disentangle in phase and energy space underlying physi- 
cal components. Assuming again that the Bridge spectrum rep- 
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Fig. 11. Scale factors of the 3 empirical spectral models as a function of phase: left panel, power-law model scale factor; middle 
panel, "narrow" bump (PI modified power-law component) scale factor; right panel, "broad" bump (Bridge modified power-law 
component) scale factor. 



resents the shape of a distinct emission component, we mod- 
elled its spectral behaviour in terms of two spectral compo- 
nents, a "modified" power-law (mpl) with an energy depen- 
dent index and a simple power-law (pi): F = F mpl + F pl = 
Q,.£H/3+7-in £ ) +a pV E~ *P ] . In this formula F denotes the pho- 
ton flux in units ph/cm 2 s MeV, while E is given in MeV. As ex- 
pected, the normalization parameter a p i for the simple power- 
law component was consistent with zero, although a weak level 
of high-energy 7-ray emission is measured up to GeV energies. 
The same approach has been followed for the emission in the 
PI interval. The resulting simple power-law component, de- 
scribing entirely the high-energy end of the spectrum, has an 
index 2.022 ± 0.014. We note that this power-law spectrum has 
to break near the boundaries of the energy window shown in 
Fig. [l(] (see also Fig. ^|). The resulting best fit values for the 
"modified" power-law components of the PI and Bridge emis- 
sions are given in Table [|. From these values we can derive the 
following positions for the maxima (in the E 2 x F flux repre- 
sentation): 14.0 ± 1.1 keV and 135 ± 15 keV for the PI and 
Bridge phase intervals, respectively. The widths of the "modi- 
fied" power-law components are specified by the FWHM val- 
ues in the 10 \og(E) domain and are about 1.81 and 2.65 for the 
PI and Bridge intervals, respectively. With this approach we 
have identified three distinctly different spectral shapes, which 
can describe the PI and Bridge spectra. 

In the next step we made an attempt to describe the mea- 
sured spectral distributions in the narrow pulse-phase intervals 
(npi) in terms of just these 3 models each with a free scaling 
parameter: F n P ; = a ■ F Brid s em P> + b ■ F pl ' m v l + c ■ E~ Ppi -& 

Interestingly, the resulting fits are very satisfactory for all 
phase intervals as shown in Fig. in which the composite 
model (dotted lines) and the individual components are su- 
perimposed on the measured spectra. The fit characteristics 
for each phase interval are shown in Table ^| The xt va l~ 
ues of the fits indicate acceptable spectral descriptions in all 
cases. The inclusion of 10-15% systematic uncertainties in the 
flux measurements, however, makes a straightforward assess- 
ment/interpretation of the \ 2 V values difficult. 



Table 4. Best fit parameter values for the "modified" power- 
law 

and power-law components of the PI phase interval and for the 
"modified" power-law of Bridge phase interval 



Component 


a 


P 


7 




(ph/cm 2 sMeV) 






PI -mpl 


(1.06 ± 0.07)E-5 


3.361 ±0.014 


0.159 ±0.003 


PI -pi 


(9.98 ± 0.69)E-5 


2.022 ±0.014 




Bridge-mpl 


(7.09 ± 0.51)E-5 


2.298 ±0.017 


0.074 ± 0.003 



In Table |^ two types of error estimates are presented for 
each scaling parameter. The first type is associated with the 
(asymmetric) statistical uncertainty in the scale parameter us- 
ing the best fit estimates for the parameters describing the shape 
of the 3 models i.e. /3 P i-m P i, TPi-mpi, /?Brid g e-m P i , 7Brid g e-m P i and 
/3pi_ p i (see Table Q). The second type is related to the system- 
atic uncertainty in the fitted scale parameter and has been deter- 
mined by varying the shape parameters of the 3 models within 
their ±lcr errors. The range over which the scale parameters 
vary is indicative for the systematic uncertainty due to uncer- 
tainties in the shape of the 3 model fit functions. 

The fit results for the 3 scale parameters are visualized 
in Fig. 
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These scale parameters have been normalized to 
the emission in the Bridge interval because of the different 
phase extents of the intervals. Fig. [IT] shows effectively the 
"light curves" of these model components: two components 
are clearly related to the emission in the two main pulses 
(the power-law component and the "narrow bump") and the 
"Bridge component" or "broad bump" extends apparently from 
the LW1 till under P2 in a triangular shape. 

The "light curves" in Fig. |ll] can also be clearly discerned 
in the 9 panels of Fig. [l(]. The upper row centered on PI shows 
how the components of PI extend into the wings; the middle 
row shows how the Bridge spectrum dominates in all three 
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Table 5. Best fit scale factors (a, b, c) for the narrow phase intervals 
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Interval a b c xl (x 2 / n dof) 

+0.024 +0.008 n nR7 +0.011 +0.008 / ro + 1.01 +0.22 

-0.029 -0.008 U - UD ' -0.009 -0.008 I z ' o>3 -1.00 -0.22 

+0.061 +0.001 i nn/ i+0.050 +0.086 /qq nQ +4.59 +2.52 

-0.061 -0.001 1 - uu *-0.047 -0.086 l yy - uy -4.93 -2.52 

+0.052 +0.017 n i n o+0. 019 +0.025 i-t n fi o+l-85 +0.39 

-0.052 -0.017 u - iu,; >-0.019 -0.025 l iu - DZ -1.86 -0.39 

+0.051 +0.033 r, n1 9+O.OI8 +0.019 / r, m +1.40 +0.07 

-0.050 -0.033 u - ulz -0.018 -0.019 V u ' ul -1.41 -0.07 

+0.048 +0.016 n na o+0.017 +0.027 / a co + 1-59 +0.35 

-0.048 -0.016 U.UOZ_ 017 _ .Q27 I D - DO -1.60 -0.35 

+0.118 +0.064 n 70C+0. 047 +0.093 ( A7 71 +4.05 +1.95 

-0.122 -0.064 u - ' ZD -0. 048 -0.093 1 ■ ' 1 -4.15 -1.95 

+0.027 +0.000 n 9 r 7 +0.016 +0.017 1 a nc + 1-53 +0.48 

-0.027 -0.000 u - zo '-0.017 -0.017 V D - yD -1.53 -0.48 



LW1 


0.130 


PI 


0.000 


TW1 


0.696 


Bridge 


0.972 


LW2 


0.849 


P2 


1.816 


TW2 


0.000 



E-6 


0.94 (69.83/(77 


-3)) 


E-6 


0.66 (53.80/(85 


-3)) 


E-6 


0.46 (38.07/(85 


-3)) 


E-6 


0.76 (62.48/(85 


-3)) 


E-6 


0.59 (48.22/(85 


-3)) 


E-6 


0.48 (39.61/(85 


-3)) 


E-6 


0.83 (61.10/(77 


-3)) 



intervals and how the power-law component and the "narrow 
bump" contributions are symmetrically distributed on either 
side of the Bridge interval. Finally, the lowest row centered on 
P2 looks very much like the upper row, but the Bridge spectral 
component reaches a maximum value in P2. We have appar- 
ently succeeded in identifying likely genuine underlying phys- 
ical components in phase and energy space. 

7.5. Enhanced high-energy ^ -ray emission in the 
LW2 interval 

The most apparent and significant (4.3er for energies above 300 
MeV) deviation from the composite fits in Fig. [Io| is visible 
in the LW2 s pectrum in the EGRET range above 100 MeV. 
Fierro (1995) reported for this phase interval the hardest pulsed 
7-ray spectrum (photon index 1.69±0.08) fitting CGRO Cycles 

0- III EGRET data. For our (timing) analysis we used CGRO 
Cycle 0-VI EGRET data, almost doubling the statistics for en- 
ergies above 1 GeV[] 

To verify whether the phase distributions in the Ge V energy 
range also show evidence for a separate hard spectral compo- 
nent in front of P2 we produced the pulse profiles of Fig. |l2[ the 

1- 10 GeV pulse profile superimposed on the 100-300 MeV pro- 
file for EGRET Cycle 0-VI observations. The two distributions 
are normalized on the PI phase interval. A clear increase of 
emission in the LW2 phase interval (0.25-0.32) is visible for the 
highest energies, but the enhancement seems to extend to the 
maximum in the P2 phase interval, however constituting a mi- 
nor fraction in the latter. The effect can be interpreted as a phase 
shift of the second pulse with increasing energy, but also as a 
"new" spectral component (phases ~ 0.2-0.4) with respect to 
the composite fits with three spectral components as shown in 
Fig. [lO . feamanamurthy (1994)| and ^ikenberry & Fazio (1997) 



discussed variations in the intra-peak phase separations as a 
function of energy. A further phase shift above 1 GeV is con- 
sistent with the reported trend. 

Considering the enhancement a new spectral component 
then its spectrum must be even harder than derived by 



|Fierro (1995) Fitting for the LW2 phase interval another 
power-law model with both a free normalization and a free in- 
dex on top of the model composed of the broad and narrow 
spectral components and the power-law component with a fixed 
index of 2.022, we find for the additional power-law model a 
very hard photon index of 1.44+QQ4. However, our data are 
not of sufficient quality to discriminate between an additional 
very hard power-law component, which must break somewhere 
above 10 GeV to be consistent with the non-detections of 
pulsed emission at TeV energies, or for example an additional 
high-energy "bump". This can be studied in detail with the 
next generation high-energy 7-ray telescopes AGILE and par- 



ticularly GLAST (see e.g. http://glast.gsfc.nasa.gov), which is 



3 for energies lower than 1 GeV the increase in statistics is less 
due to the energy and time dependent EGRET sensitivity caused by 
spark-chamber gas aging. 



~ 30 times more sensitive than EGRET over a much wider 
energy interval extending to 300 GeV. 

8. Summary and discussion 

In this work we derived final COMPTEL pulse profiles and 
spectra for the Crab pulsar and nebula at medium 7-ray ener- 
gies (0.75-30 MeV) using data collected over the 9 year mis- 
sion of NASA's Compton Gamma-Ray Observatory. 

Due to the high counting statistics over the total 0.75- 
30 MeV interval, we were able to show a clear morphology 
change of the pulse profile as a function of energy, provid- 
ing clear evidence for drastic spectral variations with pulsar 
phase over the COMPTEL energy window. Indications for such 
variations were found in the earlier COMPTEL analysis by 
[Much etal. (1995) . 

Using our large COMPTEL data base we derived an im- 
proved Crab nebula spectrum, which has a power-law spectral 
shape between 0.75 and 30 MeV with index 2.227 ± 0.013. 
Also our new COMPTEL spectrum for the total pulsed emis- 
sion (nebula/DC emission subtracted) can be described with a 
power-law spectral shape between 0.75 and 30 MeV with index 
2.24±0.04. If the indication for enhanced emission in the 10-15 
MeV interval is genuine, then the index becomes, 2.35 ± 0.06. 

These improved COMPTEL findings have been put in a 
much broader context by including in our analysis data from in- 
struments sensitive at the neighbouring X-ray/soft 7-ray ener- 
gies, particularly from the BeppoSAX LECS, MECS and PDS 
instruments, and at high 7-ray energies from CGRO EGRET. 
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Fig. 12. EGRET pulse profiles (60 bins) using CGRO EGRET 
Cycle 0-VI data: 1-10 GeV, dotted line; 100-300 MeV, solid 
line. The profiles are normalized on their emission in PI. 
Typical error bars are indicated for both profiles. A clear in- 
crease is visible in the LW2 phase interval (0.25-0.32) for the 
1-10 GeV energy interval with respect to the 100-300 MeV in- 
terval. 



We compiled a new spectrum of the Crab pulsed emission 
from optical wavelengths up to the high-energy 7-rays at 10 
GeV (Fig. ^|). This emission reaches a level of maximum lumi- 
nosity per decade in energy from ~ 5 keV to 50 keV. Beyond 
this maximum a gradual softening sets in reaching a plateau 
(photon power-law index of ~ 2) near ~ 30 MeV which con- 
tinues to ~ 10 GeV. Above ~ 10 GeV the spectrum must break 
rapidly in order to be consistent with the stringent TeV upper 
limits for pulsed emission. 

Phase resolved spectral analysis can provide important con- 
straints for pulsar modelling, particularly to help identifying 
different production mechanisms and sites in the pulsars mag- 
netosphere. Therefore, we derived consistently over a broad en- 
ergy range from 0.1 keV up to 10 GeV (BeppoSAX LECS, 
MECS, PDS, CGRO COMPTEL and EGRET) for seven nar- 
row phase intervals phase-resolved spectra. These spectra ex- 
hibited very different spectral shapes, most notably the spec- 
tra for the narrow Bridge and Peak 1 intervals. We could dis- 
entangle the pulsed emission in energy and phase space, ex- 
ploiting the vastly different spectral shapes, particularly over 
the COMPTEL energy window, by making empirical fits, and 
found that the pulsed emission can be described with 3 dis- 
tinctly different spectral components: 

-1- a power-law emission component from ~ 1 keV to ~ 5 

GeV, photon index 2.022 ± 0.014, which is present in the 

phase intervals of the two pulses. 
-2- a curved spectral component required to describe soft (< 

100 keV) excess emission present in the same pulse-phase 

intervals. 



-3- a broad curved spectral component reflecting the bridge 
emission from 0.1 keV to ~ 10 MeV. This broad spectral 
component extends in phase over the full pulse profile in an 
approximately triangular shape, peaking under the second 
pulse. 

Furthermore, in addition to the 3 spectral components the 
Leading Wing 2 (LW2) phase interval exhibited a very hard 
spectral component, most notably at GeV energies, which 
likely extends over the broader phase interval ~ 0.2-0.4. 

In a somewhat different approach and using only 
BeppoSAX data between 0.1 and 300 keV, Massaro et al. 
(2000) identified recently two components: the first is the com- 
bination of the two components as described above under -1- 
and -2- and the second corresponds to the one described under 
-3-. In their narrower energy window differences in spectral 
shapes can be well approximated by variations in power-law 
index. 

Since the discovery of 7-ray emission from radio pulsars 
in the early seventies, the most popular and competing mod- 
els attempting to explain the high-energy radiation from highly 
magnetized rotating neutron stars can be divided in two distinct 
catagories: the so-called Polar Cap (PC) models and Outer Gap 
(OG) models. 

Detailed information on the PC models can be found in e.g.: 
Daugherty & Harding ( 1982,1994,1996), Stumer & Dermer 
(1994) and most recently |Zhang & Harding (2000)] PC models 
have problems in explaining the overall measured Crab pulsar 
characteristics. Most notibly, the large angles of ~ 60° esti- 
mated for both the magnetic inclination a and the viewing an- 
gle £, the angle between the spin axis and the observer's line of 
sight, from radio and optical/UV observations cannot be rec- 
onciled (see for recent publications e.g. Graham-Smith et al. 
1996, Moffett & Hankins 1999 and Everett & Weisberg 2001). 

OG models have no difficulties with the large a and 
£ angles estimated for the Crab, as is clearly shown by 
|Chiang & Romani (1994)1 In these models the acceleration of 
charged particles and production of high-energy radiation takes 
place in charge depleted gaps between the null-charge surface, 
defined by ft • B = with B the local magnetic field, and the 
light cylinder (with radius Ri c = c/f2) above the last closed 
field lines. For early papers and lat er refinements see: Cheng 
et al. (1986a,b), |Ho (1989)1 |Chiang & Romani (1994) , 
fcomani & Yadigaroglti (1995)1 jRomani (1996)| and 
Yadigaroglu (1997). Recently, |Cheng et al. (2000)| pre- 
sented a three-dimensional outer gap model building on the 
work of Romani and co-workers. The emission patterns from 
these outer gap models resemble fan beams, and double 
peak profiles with (strong) bridge emission can commonly 
be generated f or the cases tha t emission is seen from only 
one pole, e.g. Romani (1996) as well as from both poles, 
|Cheng et al. (2000)1 

Early attempts to model the Crab Total Pulse spectrum in an 
OG scenario using various radiative processes like the curva- 
ture, synchrotron and inverse Compton radiation mechanisms, 
assumed to play a key role in th e outer gap physics, were 
made by Cheng et al. (1986b) and |Ho(1989)| The latter em- 
ployed a self-consistent iterative procedure with one varying 
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parameter, the ratio of gap height and curvature radius of the 
field lines, and it is interesting to note that the calculated high- 
energy spectrum bears reasonably good overall simularity with 
the observed one (cf. Fig. 4 of Ho ( 1 989^ with Fig. ^ of this pa- 
per). Ulmer et al. (1995) compared the outer gap model of Ho 
with an early CGRO spectrum of the Crab (combining OSSE, 
COMPTEL and EGRET spectra) and found also good over- 
all agreement. Chiang & Romani (1994) made refinements to 
the above calculations and attempted to model the Crab Total 
Pulse spectrum and the spectral variation with phase, which 
they considered to be a clear mapping of location in the mag- 
netosphere to pulse-phase. They divided the outer gap in dif- 
ferent sub-zones taking into account the transport of radiation 
and particles from sub-zone to sub-zone. Convergence to a self- 
consistent solution, however, resulted in spectra significantly 
lacking photon flux below several GeV. Romani (1996) de- 
scribed a revised picture of gap closure and radiation physics in 
the outer magnetosphere to overcome difficulties in the schema 
of Cheng et al. (1986b), and also addressed spectral variations 
with pulsar phase from the optical to the high-energy 7-ray 
spectrum. For the Crab pulsar, he made some qualitative state- 
ments on the expected spectral properties. He expects a signif- 
icant contribution of synchrotron photons to the high-energy 
7-ray flux. This could probably explain the observed underly- 
ing power-law component from soft X-rays up to high-energy 
7-rays, although the observed photon index of ~ 2 is consider- 
ably softer than the expected value of ~ 1. 

Building on the work of Romani and co-workers, 



Cheng et al. (2000) also used a three-dimensional pulsar mag- 
netosphere to study the geometry of outer magnetospheric gap 
accelerators. However, the physics of both models is strikingly 
different. For the single outer gap model (Chiang & Romani 
1994), the emission comes from the outward direction in an 
outer gap above one pole; the emission regions for the two 
peaks of the pulse profile are those close to the null-charge sur- 
fac e and to the light cy linder radius, respectively. In the model 
of Cheng et al. (2000) photon emission consists of emission 
outward and inward from regions in outer gaps above both 
poles, the gaps being limited along the azimuthal direction by 
pair production of inward-flowing photons from the outer 
gap. It is shown that both models can pr oduce the same (Crab- 
like) pulse profiles. |Cheng et al. (2000) also calculated phase- 
resolved spectra of the Crab pulsar. They determined the lo- 
cations of the emission regions in the outer gaps in the open 
field line zone of the Crab magnetosphere as a function of 
pulse-phase. It can be seen in their Fig. 9 that high-energy 
emission from the PI interval is produced high in the mag- 
netosphere (0.8 < r/Ri c < 1.0) where curvature radiation 
dominates, resulting in a spectrum which extends to the GeV 
regime. In the interval between the pulses (TW1 and Bridge in- 
terval) high-energy radiation is predominantly produced deep 
in the magnetosphere where a soft synchrotron component is 
expected to dominate, roughly in accordance with the observa- 
tions. Moving towards the P2 interval, emissions from regions 
high and low in the magnetosphere contribute, resulting in a 
overall spectrum composed of a hard curvature component and 
a much softer synchrotron component. In the P2 interval high- 
energy radiation is coming, in essence, from emitting regions 



extending from ~ 0.2i?/ c to ~ 1.0i?; c which gives rise to a 
hard spectral component extending into the GeV domain and a 
soft component. Crossing the bridge interval moving from PI 
to P2 a gradual decrease from ~ 0.6i?; c to ~ 0.2i?; c is seen 
for the lower bound of the emission region in the outer gap. 
According to Eq. 33 of Cheng et al. (2000) the dominating syn- 
chrotron emission from these regions deep in the outer gap in 
the pulsar's magnetosphere becomes increasingly intense mov- 
ing towards P2, because the magnetic field strength becomes 
stronger deeper in the magnetosphere. This offers an expla- 
nation for the observed phase dependence of the bridge/broad 
bump spectral component, shown in Fig. [n]. Beyond P2 (and 
before PI) this soft synchrotron component should be absent 
which is in agreement with the observations. Thus the model 
proposed by Cheng et al. (2000) seems to provide a viable and 
promising theoretical description of the physics responsible for 
the production of the Crab high-energy radiation with char- 
acteristics as shown in Figs. ^ [To] and |ll|. A direct quanti- 
tative confrontation of this model with our observed 0. 1 keV 
- 10 Ge V phase-resolved sp ectra is therefore strongly recom- 
mended. |Cheng et al. (2000*) compared in their paper the model 
calculations with phase-resolved Crab spectra from EGRET 
(> 30 MeV, Fierro et al. 1998). By mistake, EGRET spectra 
from the phase-resolved spatial analysis were used, which also 
comprise the underlying nebula component. This explains the 
large discrepancies between the model calculations and the ob- 
served spectra for energies between 30 and 100 MeV. Cheng et 
al. (2000) also show a broad-band (0.1 keV - 10 GeV) model 
spectrum for the phase-averaged Crab pulsar spectrum, which 
can be compared with our spectrum in Fig. ^. The overall shape 
of the model spectrum follows the observed spectral character- 
istics well, although for energies below ^100 keV the model 
underestimates the observed X-ray fluxes, i.e. it seems that the 
component which we empirically described as a narrow spec- 
tral "bump" peaking around 20 keV in an E 2 x F representa- 
tion, is not accounted for. This can be best studied in the phase- 
resolved analysis. 

Future observations of the Crab pulsar by high-energy 
missions are important. In particular the spectral characteris- 
tics in the 300-1000 keV interval must be determined much 
more accurately. Here, data from both IBIS and SPIE aboard 
INTEGRAL will contribute significantly. At high 7-ray ener- 
gies data from the AGILE and GLAST missions can provide 
sufficient statistical precision to study the peculiar spectral be- 
haviour in the LW2 phase interval in much more detail for ener- 
gies above 1 GeV. Moreover, these instruments can, for the first 
time, study the pulsed emission for energies above 10 GeV (up- 
per bound of the sensitivity window is 50 GeV and 300 GeV 
for AGILE and GLAST, respectively). At medium 7-ray en- 
ergies there are plans for a mission with an advanced, more 
sensitive Compton telescope (called MEGA 0.5-50 MeV; see 



Kanbach 2001 ), but no approval for such a mission exists at this 



time. It is just in this energy range, where interesting spectral 
transitions occur and where we have indications for enhanced 
pulsed emission in the 10-15 MeV range. Future space borne 
Compton telescopes having 10-100 times better sensitivity than 
CGRO COMPTEL are required to allow further progress. 
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Appendix A: Uncertainties in the absolute flux 
measurements at hard X-rays/soft 
7-rays 

In this work we compiled Crab spectra from soft X-rays up to 
high-energy gamma rays in order to better determine its spec- 
tral characteristics. Meanwhile, the Crab is used by many in- 
struments as an in-flight calibration source. However, how well 
do we know the genuine Crab spectrum? At lower X-ray en- 
ergies the measurement of the spectrum is coupled to an esti- 
mate of the Nh value, at harder X-rays up to the gamma-ray 
regime we have direct measurements of the Crab spectrum, 
but how realistic are the estimates of the accuracy of the pre- 
launch instrument calibrations? Is there any for which we can 
trust the claimed accuracies most? In this Appendix we show 
the present status of our knowledge of the Crab spectrum. The 
still large (see below) systematic differences between different 
instruments should be kept in mind, not only when drawing 
conclusions on Crab results, but for any other source in high- 
energy astrophysics. 

we derived a N H value of 3.61(2) x 10 21 cm~ 2 . 



In Sect. 7.1 



This estimate is significantly larger than the Nh value of 
3.23(2) x 10 21 cm" 2 derived recently by Massaro et al. (2000) 
analysing a much larger BeppoSAX Crab database. The ap- 
parent discrepancy can be explained by their use of an older 
version of the LECS response description in combination with 
fitting LECS 0.1-4 keV data only, thus constraining to a lesser 
extent the photon index. We verified this by reproducing our 
values for Nh, the normalization and the photon index for this 
Off Pulse emission repeating our analysis for the database used 
by Massaro et al. (2000). From this comparison it is clear that 
the uncertainty in Nh is dominated by systematic uncertainties 
in the response characteristics, of the LECS in particular, rather 
than by statistical ones. 

Our value for Nh of 3.61(2) x 10 21 cm~ 2 is consis- 
tent with th e X-ray based value of 3.45(42 ) x 10 21 cm~ 2 
obtained by |Schattenburg & Canizares (1986)| using the Focal 
Plane Crystal Spectrometer on the Einstein observatory and 
with a radio based estimate of ~ 3.65 x 10 21 cm -2 
( [Dickey & Lockman 1990} ). Recently, a study using XMM 
Newton EPIC MOS data yields a hydrogen column density of 
3.45(2) x 10 21 cm -2 for an oxygen-ir on depleted abundance 



of 0.63(1) solar ( |Willingale et al. 2001 ). Assuming solar abun 



dances the N H value lowers to 3.28(2) x 10 21 cm 2 . These 
estimates indicate that the genuine value of Nh will very prob- 
ably lie in the range (3.3 - 3.6) x 10 21 cm" 2 . 

In the fit to the data of the 4 BeppoSAX NFI instruments 
we derived for the sensitivity normalization scale factors, rela- 
tive to the MECS (factor set to 1): LECS 0.93, HPGSPC 1.01 
and PDS 0.87. This means that the LECS and PDS calibrations 
of their overall sensitivities deviate from that of the MECS by 
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Fig. A.l. The Crab nebula spectrum in the 1 keV - 10 MeV en- 
ergy interval. Flux measurements from BeppoSAX LECS (1- 
4 keV), MECS (1.6-10 keV), HPGSPC (10-32 keV) and PDS 
(15-300 keV), GRIS (20-500 keV) and COMPTEL (0.75-10 
MeV) are shown. The dotted lines show the best power-law 
fits for the combined BeppoSAX instruments and for GRIS. 
Note the apparent systematic deviation from the power-law 
fit in the LECS spectral data. The shaded band indicates the 
±1(7 uncertainty interval around the optimum power-law fit for 
COMPTEL (0.75-30 MeV). Clear discrete jumps are visible 
between the various BeppoSAX instruments and GRIS reflect- 
ing uncertainties in absolute sensitivity. 



7% and 13%, respectively. This is clearly visible in Fig. \A.l[ 
which presents the nebula spectrum in an E 2 x F represen- 
tation between 1 keV and 10 MeV as measured with the 4 
BeppoSAX NFI instruments and COMPTEL, in which no nor- 
malization correction factors have been applied. Also shown 
is the Crab nebula spectrum in the 0.02-1 MeV energy range 
as measured by the balloon borne GRIS (Ge detectors; Bartlett 
et al. 1994a). The combined BeppoSAX spectra as well as the 
GRIS spectrum are best fitted with a power-law spectral shape 
with a consistent slope of ~ 2.14 over the 1-700 keV inter- 
val. However, the normalization factor for GRIS, relative to 
the MECS is even as low as 0.78, to be compared with the 
estimated GRIS systematic uncertainties of +12% and —6% 
by Bartlett (1994b). Allowing spectral curvature in the multi- 
instrument BeppoSAX nebula fit by introducing an energy de- 
pendent power-law index does not improve the fit significantly. 
The same is true for the GRIS spectrum. However, some grad- 
ual softening above a few 100 keV is required to connect to 
the softer spectru m m easured by COMPTEL at energies above 
1 MeV (cf. Fig. Al; power-law photon index in the 0.75-30 
MeV interval is 2.227 ± 0.013). The statistical uncertainties in 
the above quoted normalization correction factors are typically 
better than 1%. Therefore, it is obvious that the differences in 
absolute normalizations are systematic, and it is discouraging 
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to note that between the two instruments for which the calibra- 
tions were expected to be most accurate (MECS and GRIS), 
the discrepancy appears to be largest. We cannot decide un- 
ambigously on this controversy. Therefore, we feel that in the 
presentation and analysis of the above spectral data in com- 
bined broad-band spectra, this problem should not be hidden 
by making arbitrary choices on normalization. 
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